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Abstract 

A kinetic theory of vacuum particle creation under the action of an 
inertia! mechanism is constructed within a nonpertrubative dynamical 
approach. At the semi-phenomenological level, the inertial mechanism 
corresponds to quantum field theory with a time-dependent mass. At the 
microscopic level, such a dependence may be caused by different reasons: 
The non-stationary Higgs mechanism, the influence of a mean fleld or 
condensate, the presence of the conformal multiplier in the scalar-tensor 
gravitation theory etc. In what follows, a kinetic theory in the colli- 
sionless approximation is developed for scalar, spinor and massive vector 
fields in the framework of the oscillator representation, which is an effec- 
tive tool for transition to the quasiparticle description and for derivation 
of non-Markovian kinetic equations. Properties of these equations and 
relevant observables (particle number and energy densities, pressure) are 
studied. The developed theory is applied here to describe the vacuum 
matter creation in conformal cosmological models and discuss the prob- 
lem of the observed number density of photons in the cosmic microwave 
background radiation. As other example, the self-consistent evolution 
of scalar fields with non-monotonic self-interaction potentials (the W - 
potential and Witten-Di Vecchia - Veneziano model) is considered. In 
particular, conditions for appearance of tachyonic modes and a problem 
of the relevant definition of a vacuum state are considered. 
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1 Introduction 



The present work is devoted to the construction of a kinetic theory of vacuum 
creation of particles with time-dependent masses. For brevity, this mechanism 
will be referred to as incrtial one. Microscopic foundations of a mass change 
may be different. The Higgs mechanism leads to the most popular models of 
such a class, when the corresponding mean fields are time-dependent. General 
quantum field models with nonpolynomial interactions may also be considered, 
where the separation of nonstationary mean fields results in a time-dependent 
mass |T]. A well-known example of this kind is the Witten - Di Vecchia - 
Veneziano model [5J [3] in the framework of which the mean-field concept was 
analyzed in [J . The Nambu- Jona-Lasinio [5J and a [5] models are other exam- 
ples, where the meson masses are defined by evolution of a quark condensate to 
be described at the hydrodynamic 7J or kinetic 8J level. The particle mass may 
depend on many-particle interactions in hot and dense non-stationary matter 
[HI Uni E] • A general basis for a rather slowly- varying time dependence of the 
effective mass can be obtained within the Green function method [T^l [13] . The 
field dependence of the mass is a general factor determining the time evolution 
in all these cases (F-class models). The conformal invariance of the scalar- 
tensor gravitational theory provides a time dependence of the particle mass by 
means of the conformal multiplier [HI [TH [161 [E] • The mass can be changed 
also due to the parametrization stipulated by additional space dimensions fl8j . 
Such theories should be referred to as the other class (C-class). In the F-class 
theories, the vacuum particle creation admits a well-known interpretation based 
on the simplified vacuum tunnelling model in an external field [l9l [20l [2T| . A 
similar interpretation of the C-class models is difficult. At the phenomenolog- 
ical level, however, both classes have a uniform mathematical description as 
it will be shown in Sects. [21 [31 and [H respectively, for the scalar, fermion and 
massive vector boson in quantum field theories (QFT's). 

The first consideration of the vacuum creation of particles with the vari- 
able mass was proposed apparently in j22| as a possible variant for describing 
a quantum system response to the time variation of system parameters [23]. 
Using the Bogoliubov transformation method, residual momentum distribu- 
tions for fermions, and pair correlators were found for the cases of step-like and 
smooth variations of the fermion mass (Sect. 13. 4p . 

In the present work, the kinetic theory will be based on the oscillator repre- 
sentation (OR) [Ml, , which is the most economical method for a nonpertru- 
bative description (as compared to the Bogoliubov method of canonical trans- 
formation [5S] or other accurate approaches to the problem [TS] [571 HH] ) of the 
vacuum particle creation under the action of time-dependent strong fields. This 
approach leads directly to the quasiparticle representation (QPR) with diago- 
nal operator forms in the momentum space for the set of dynamical variables. 
It allows one to get easily the Heisenberg-type equations of motion for creation 
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and annihilation operators. An important feature of the time-dependent Fock 
representation is the necessary consistence of commutation (anti-commutation) 
relations with the equations of motion. Otherwise, this circumstance can bring 
to the non-canonical quantization rules (an example will be considered in Sect. 
SH). 

In terms of the OR it is possible to immediately derive the corresponding 
kinetic equations (KE) by the well known method . Some particular results 
are published in [TJ[3D] . The kinetic theory for scalar, spinor, and massive vector 
fields is constructed in Sects. [H [3] and [H respectively. The main attention is 
paid to the particle creation in conformal cosmological models [311 |3H [33] 
(Sec. [S]). It is shown that the choice of the equation of state (EoS) of the 
Universe allows one to obtain, in principle, the observed number density of 
matter participants and photons and, possibly, dark matter. The basic problem 
here is the description of vacuum particle creation which should be consistent 
with EoS but it is beyond the present article. 

Finally, in the Sect. [5] other class of scalar QFT systems is considered 
with non-monotonic self-interaction potentials to apply the decomposition of 
the field amplitude into the quasiclassical space-homogeneous time-dependent 
background field and the fluctuation part. In this case, the particle mass is 
defined by intensity of a quasiclassical field. As an example, self-interaction 
potentials of the simplest polynomial type and that for a nontrivial case (the 
pseudoscalar sector of the Witten - Di Vecchia - Veneziano model) are ana- 
lyzed. It is shown, that the relevant definition of vacuum states allows one 
to avoid of the tachyonic mode beginning. The main purpose of this review 
is to summarize all known relevant results on the inertial mechanism of the 
vacuum particle production and to call attention to unsolved problems which 
are shortly Hsted in Sect. [7] 

We use the metric 5^'' = diag(l, —1, —1, —1) and natural units h — c — 1. 

2 Scalar field 

2.1 Oscillator and quasiparticle representations 

Let us start our consideration with the simplest case of the real scalar field 
with the time-dependent mass m(t), whose equation of motion is 

[df,d'' + m\tMx)^0. (1) 

The corresponding Lagrange function is given as 

L=^d^^d^^-^m^it)^'. (2) 

In the considered case, the system is space-homogeneous and nonstation- 
ary. Therefore, the transition to the Fock space can be realized on the basis 
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functions y(x) ~ exp(±jpx), and creation and annihilation operators become 
time-dependent. The assumption about the space homogeneity aUows one to 
look for solution of Eq. ([1]) in discrete momentum space in the following form: 

P 

where V = and pi — (27r/L)ni at that the integers rii, (i = 1,2,3) run 
from — oo to +oo. The thermodynamic limit can be covered in the resulting 
equations. Then the oscillator-type equation of motion follows from Eqs. ((T]) 
and decomposition ^ as 

(^W(p,i)+c.2(p,tV(±)(p,i)=0 (4) 

with 

w2(p,i)=m2(t)+p2 . (5) 

The symbols (±) correspond to the positive and negative frequency solutions 
of Eq. (m defined by its free asymptotics in the infinite past (future) [55] , 

^W(p,f^Tcx3) ~e±'^-^*, (6) 



where Wip = y + p^ arc defined by asymptotics of the mass 

= lim m[t). (7) 

>=poo 

The asymptotics ([6]) corresponds to the in(out)-states and is necessary for def- 
inition of in(out)-vacuum. This requirement, however, can be broken in cos- 
mology [34]. We suppose here that such asymptotics exists and the relevant 
vacuum states will be denoted by |0) without indices "in" or "out" , that is evi- 
dent from the context. In the considered class of problems, the classification of 
states in the frequency sign turns out to be impossible for an arbitrary time mo- 
ment. According to a general analysis [35', this leads to instability of a vacuum 
state during the action period of the external fields and to the vacuum particle 
creation. In this case, it is possible to consider quasiparticle excitations during 
the system evolution {m{t) ^ 0) and describe their creation and annihilation 
in the vacuum state |0). When the external field action is completed, residual 
particles of some finite density remain in the out-state. However, it is necessary 
to emphasize that these particles are defined in respect of the in- vacuum state 
[26]. 

In the general case, S-matrix does not exist in the considered formalism. Its 
role in description of the system with unstable vacuum is performed by other 
mathematical objects: The operator of the canonical Bogoliubov transforma- 
tion [26] , distribution functions in the kinetic approach [24l |29l jSO] |36l jST] [38] 
or some set of correlation functions |39l |40| etc. 



4 



The conception of "quasiparticle" plays the central role in the QFT with 
strong time-dependent quasiclassical external fields [26l [35] . Under the consid- 
ered conditions, this approach is the shortest realization of the quasiparticle 
concept within the standard Fock representation of the QFT, where the external 
field can be taken into account nonperturbatively. Thus, the QPR corresponds 
to a possibility of writing down the set of commutative operators of physical 
(observable) quantities (the complete QPR) in the diagonal form in an arbi- 
trary time. It is naturally related to the question whether operators have the 
quadratic form in the Fock representation. Hence, the interaction between the 
field constituents and self- interaction is not taken into account. It corresponds 
to a nondissipative approximation in the kinetic theory [41j . An alternative 
definition of the quasiparticle was given in [25j for constrained systems. 

The transition to the QPR can be realized in different ways. The traditional 
method is based on the time-dependent canonical Bogoliubov transformation 
p6] . The alternative approach uses the oscillator (" holomorphic" ) representa- 
tion, which leads directly to the QRP [25 • In the considered case the transition 
to the OR is made by substituting m± m(t) into the dispersion law for the 
free field and postulating a following decompositions: 

V2V ^ V^(p,i) I J 

7t{x) = --L= ^ v/^J(F^ {a(p, t) e^P- - aHp, t) g-^P^} , (8) 
V2V p 

where 7r(x) is the generalized momentum, a'^{p,t) and a{p,t) are the creation 
and annihilation operators of particles with the momentum p at the time mo- 
ment t. The in-vacuum state is defined as 

a(p,t ^ -oo)|0 >= 0, <0|0>=1. (9) 

The canonical commutation relation 

Mx),7T{x%^,,=tS{^-^') (10) 

together with the decomposition ([5]) provides the standard commutation rela- 
tion for time-dependent creation and annihilation operators 

[a{p,t),a\p',t)]^Spp, . (11) 

The substitution of the decompositions ^ into the Hamiltonian 

H{t) = 1 fd''^ + [V<p(x)]2 + m^{t)ip^{x)} (12) 
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leads immediately to a diagonal form which corresponds to the QPR 



H{t) = ^c^(p,t) |a^(p,i)a(p,i) + i| 



(13) 



In the considered case the vacuum energy of zero oscillations (" Zitterbewe- 
gung" ) depends on time. 

Equations of motion for the operators a,a^ can be obtained now from the 
minimal action principle [24] or from the Hamiltonian equations 

• . SH ^ 2,, 

(fi — — TT, TT = --— = /\(f - m (t)(p . (14) 

OTT 0(p 

Here and below we use the notation (') = d/dt{ ). Combining Eqs. ^ and 
(O, we get 

a{p, t) = ^A(p, t) a^(-p, t) - iw(p, t) a(p, t), 
at(p, t) = iA(p, t) a(-p, t) + iLo{p, t) at(p, i), (15) 

where 

^(P'*) = ^f:7)-^j^ 

is the factor defining the mixing of states with positive and negative energies. 
This equation obviously is consistent with the commutation relations (jlip . 

Equations of motion (fT5)) can be rewritten as the Heisenberg-type equation, 
e.g., 

1 

2 

In the instantaneous QPR these equations serve as a basis for a nonperturba- 
tive derivation of the KE describing scalar particle creation and annihilation 
processes within the inertial mechanism. 

2.2 Kinetic equation 

The key object of the kinetic theory is the quasiparticle distribution function 
which for the space-homogeneous case is 

/(p,t) = (0|at(p,t)a(p,t)|0), (18) 

where |0) = |0)i„ is the initial (t — oo) vacuum state. Differentiating the 
distribution function (|18p with respect to time and using (|15|) we get 

/(p,t) = A(p,t)Re{/W(p,t)}. (19) 



a(p,t) = ^A(p,t)at(_p,t) + z[i7(i),a(p,t)] . (17) 
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Here the auxiliary correlation function is introduced 



/W(p,i) = (0|at(p,t)at(_p,f)|0). (20) 

This function provides a coherent connection between the states with positive 
and negative energies (so-called entangled states [?5]). 

The equation of motion for /^^-'(p,t) can be obtained by analogy with 
equation ([1^)) . We present it here in the integral form as 

t 

/W(p,t) = \j dt'A(p,i')[l + 2/(p,t')]e2^«(P^*-*'), (21) 

to 

where the initial condition (p, to) — was used. This condition corresponds 
to the initial condition for the distribution function /(p, to) = and is a direct 
consequence of the definition ^ . Eventually, the dynamical phase in Eq. (PT|) 
is equal to 

t 

0(p;i,i') = j d™(p,r). (22) 



The substitution of Eas. ((2T|) in Eq. p9|) leads us to the resulting KE written 
in the thermodynamical limit V ^ oo aX, the fixed particle density 

t 

/(p, t) = i A(p, t) j dt' A{p', t) [1 + 2/(p, t')] cos [20(p; t, t')] . (23) 

to 

The source term in the r.h.s. of Eq. ([23]) describes a variation of the particle 
number with the given momentum due to vacuum creation and annihilation 
processes for the inertial mechanism, A(p, i) is defined by Eq. (fT6|) . The non- 
Markovian KE (|23|) has the structure as that for the Schwinger mechanism of 
pair creation in an electric field |29j . This equation was investigated in detail 
for description of the pre-equilibrium evolution of quark-gluon plasma created 
in collisions of ultrarelativistic heavy ions [36l f37] • The case of the scalar QED 
was considered in 24J for an electric field of arbitrary polarization. 

As it follows from Eq . (|16p . in the framework of the inertial mechanism the 
particle production rate is defined by the rate of the mass change 

^ m{t) d{mot) ' ^ ^ 

where ttiq is the characteristic mass to fix the time scale {e.g., mo = mJ). 
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Let us note that the KE (^5)) is valid under two basic assumptions: a) there 
are no particles (or antiparticles) in the in-state; b) a collisionless approxima- 
tion is applicable {i.e. the corresponding dissipative processes are not taken 
into consideration). 

In the low-density approximation /(p, t) 1 the KE results in the 
following solution [35] : 



dt'A(p,t')exp [2i9{p;t,t')] 



> 0. (25) 



The KE (221) can be transformed to linear equations of the non-Hamiltonian 
dynamical system with zero initial conditions 

/ ^Au, M = A(l + 2/) - 2ujv, V = 2oju, (26) 

which is convenient for numerical analysis. This equation system has the first 
integral 

{l + 2fy -u^ -v'^ = 1, (27) 

according to which the phase trajectories are located on the two-cavity hyper- 
boloid with top coordinates f = u = v = (physical branch) and f — —l,u = 
V = (nonphysical one). If the function / is excluded from Eqs. (|26p . we 
obtain the non-linear two-dimensional dynamical system with 



u = A\/ 1 + u'^ + v'^ — 2Ltjv, 

V = 2uju. (28) 

The functions /{p^t) and u{p,t) have a certain physical meaning (the last 
function describes vacuum polarization effects, see Sect. 12. 3p below and are 
invariants with respect to the time inversion t —t while the auxiliary function 
v{p,t) and factor pB)) change their signs. Thus, the KE (^5)) is invariant at 
the time inversion. 

The presented formalism of vacuum particle creation is specific for kinetic 
theory and allows one natural generalization to the case of interacting fields 
that leads to introduction of corresponding collision integral of a non-Markovian 
type [12] ■ This approach is close to the modern method expounded in the 
book [26] where the time-dependent Bogoliubov transformation is used. The 
same method was used in pioneer works [43l|44l|45] (see also [34]). Some modi- 
fication of the formalism [26] (the r, 9 representation) was developed in [25] and 
then used widely (e.g., in [32j and references cited there). The correspondence 
between the r, 9 representation and our approach (as well as that used in the 
book [2S]) may be easily established. 
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2.3 Observable and regularization 



The KE describes the vacuum quasiparticle excitations rising at an exter- 
nal force {Th{t) ^ 0, in the considered case). When this action is switched off, 
there is still some remaining density of real (residual) particles and antiparti- 
cles. In the absence of any interaction between the system constituents, the 
real particles are "on-shell" ones and have the free-particle dispersion law uj± 
with the mass m± while quasiparticles are "off-shell" with the dispersion 
law ([5]). Within the Green function method [T^, one can say that the time- 
dependent dispersion law like (O corresponds to the t-parametric mass shell 
surface of slowly time-dependent m{t), i.e. m{T ± t/2) « m{T), where T and 
T correspond to slow and fast time scales. This case is not of interest for the 
considered problem. Thus, the dispersion law ^ does not belong to the mass 
shell surface. In the general case, the on-shell condition 



(mo = TO±) is not connected directly with the condition of efficiency of the 
vacuum particle creation, ^{t) < 1, where £^{t) is defined by Ea. ([M)l . On the 
contrary, the presence of high frequencies in the function m(t) is necessary 
for vacuum creation and does not contradict the on-shell condition In 
principle, the KEs of such a type are designed for the description of evolution 
of both real particles and quasiparticles. In particular, the distribution function 
of residual particles is /out(p) = hm f{p,t). This simple formula for fout{p) 

t — >oo 

follows from Eq. (j25p in the low-density approximation. However, the presence 
of the fast oscillated multiplier in the source term in the r.h.s. of the KE 
((23|) leads to a large amount of numerical calculations which make impossible 
the study of the system evolution for rather large times after the switching on 
external forces. The corresponding large scaling methods of calculations based 
on the KE (|23p have not been worked out at present. Some properties of a 
residual particle-antiparticle plasma due to a limited pulse of the external field 
action can be estmated by the imaginary time method [46j . 

The distribution function is the key quantity of the system. The density 
of observable variables is some integral in the momentum space containing the 
distribution function and auxiliary functions u(p,t), which describe the effects 
of vacuum polarization. The simplest variable of such a type is the density of 
quasiparticles. In the thermodynamic limit L ^ oo we have 



m{t) — mo 



< 1 



(29) 




(30) 
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where [dp] = (27r) ^d^p. To proceed to the thcrmodynamical limit, the rule 

^E-/h (31) 



is used here and below. 

Other important characteristics of the system are the energy density e 
and pressure P, which can be obtained as the average value of the energy- 
momentum tensor corresponding to the Lagrangian density 



As the result we have ESI 



3P 



(OlToolO) = J[dp]L^f, 

[dp] — (/ + -u) +UJU 

UJ z 



(32) 

(33) 
(34) 



The last two terms in integrand 
polarization. 

Finally, the entropy density can be introduced 



represent the contribution of vacuum 



S{t) 



H[/ln/-(l + /)ln(l + /)]. 



(35) 



It is not conserved {S{t) ^ 0) even in the considered non-dissipative approx- 
imation because the system is open (the mass change is defined by external 
causes) . 

A direct proof of the convergence of integrals pO)) . p3|) . (p4)) . ((35)) is com- 
plicated because of the absence of an explicit form for functions /(p,t) and 
M(p,t). Therefore, one usually uses the method of asymptotic expansions in 
power series of the inverse momentum p~^ (N-wave regularization technique) 
[17] (another approach rests on the WKB approximation [IH]). Our present 
consideration is based on the explicit asymptotic solutions of the system (1261) 
for IpI 3> m. The integral (|30p is assumed to be convergent at any time moment. 
Then the function f{p,t) should decrease at p ^ oo and hence /(p,t) ^ 1 
in this region. This inequality corresponds to the low-density approximation 
j5l), where the KE solution can be written in the explicit form as 



^ 4p4 



dt'm{t')m{t') exp [lip (t - t')] 



(36) 



p = IpI ^ oo because A°°(p,t) = ■m{t)m{t) / p'^ in accordance with Eq. ([T| 
These solutions are consistent with the integral of motion (P7| . Thus, indeed 
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asymptotic solutions are some quickly oscillating functions (this fact was first 
noted in [S^ for the case of massive vector bosons, see Sect. 14. 3p . Such be- 
havior matches with the quasiparticle interpretation of vacuum excitations by 
the inertial mechanism. The real (observed) particles are the result of the 
evolution by the moment when m = and the out-vacuum state is realized. 
The asymptotics p6p may be influenced by other (non-inertial) mechanisms 
of vacuum particle creation {e.g., in the case of harmonic "laser" electric field 

The asymptotics of the integral can be obtained by the stationary 
phase method [S5] 

t 

dt' mit')m{t')e^'P^'-''^ = m{t)m{t) ^ ^^^^ 

ip 
to 

if m{to) = 0. Using Eqs. and (|^ we get the leading contributions 



m 



{t)m{t) 



2p3 

TO^(t) -I- m(t)fh{t) 



(38) 



1 

where the upper indices show the inverse momentum degree for the correspond- 
ing leading terms (we are not interested in the asymptote of the function u(p, t), 
which plays some auxiliary role only) . Relations ([55]) are identical to the results 
of application of the N-wave regularization method to Eqs. (P5|) [TT] . 

Now one can conclude that the integral (|33p is convergent but the last 
integral term in Eq. (j34p needs a regularization. The regularizing procedure of 
the Pauli-Villars type is based on the subtraction of appropriate counterterms 
in integrals ([30| . ([33ll . ([34|l . ([35| . These counterterms can be obtained by the 
substitution p^ ^ p^ ^ Af ^ into the denominator of asymptotics ([38|l , 



fR = .f-fM, Ur = U-UM- (39) 

If the regularizing mass M ^ m{t) can be chosen rather large, M ^ A (A is 
"the computer cut-off parameter"), the influence of counterterms on the results 
of numerical calculations is negligible. 

The numerical investigation of the KE (p3)) and observable densities ((30|) . 
(I33l)-(l35l) will be presented in Sect. [331 



3 Fermion field 

3.1 Quasiparticle representation 

The material of this subsection is based on papers 
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Equations of motion for fcrmion fields with the variable mass are 

[z7^9,, - m{t)]ij{x) = 0, 
^{x)[i-f^1) ^ + m{t)] ^ , (40) 

where "ip = 7/;^7°. The corresponding Hamiltonian is {k — 1,2, 3) 

H{t) = i [ (fx = ( (fx 4j{-i-i^dk + m{t)}ij}. (41) 



By analogy to the scalar case, we use the following decompositions of field 
functions in the discrete momentum space: 

'^(^) = ^E E {e'P"«o(p,t)u"(p,i)+e-^P'^6t(p,i)v"(p,t)}, 

^(^) = ^E E {e-'P"4(P,t)a"(P,i) + e'P''&a(p,i)v"(p,t)}. (42) 

The OR is intended to derive equations of motion for creation and annihilation 
operators. It is based on the primary equations (|40p and free u, v - spinors 
with the substitution m m{t). Thus, the following equations for the spinors 
are postulated in the OR : 

[yp- m{t)]u{p,t) = 0, 

[7P + m(<)]v(p,<) = (43) 

with p° = uj{p,t). These definitions create the set of standard orthogonality 
conditions [56j depending on time now parameterically 

u"(p,i)u^(p,t) = ^P%<5.^, v"(p,i)v'^(p,t) = 

ut"(p,t)u'='(p,<) = vt"(-p,i)v^'(-p,<) - 5^0, 

u"(p,t)v^(p,i) = ut"(p,i)v'3(-p,t) =0. (44) 

Decompositions (H^ and relations (jH]) lead immediately to the diagonal 
form of the Hamiltonian (l4T]l 



if(i) =^w(p,t)[al(p,t)a„(p,t) + 6l(p,i)6„(p,t)] (45) 

p,a 

with interpretation of , a (and , h) as the creation and annihilation oper- 
ators of quasiparticles obeying the standard anti-commutation relations 

{a„(p, i), a)j(p', t)) = {6,(p, t),bUp', t)} = 5pp,S^p. (46) 
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We are not interested in the subsequent diagonalization of the spin operator 
and such QPR can be named the incomplete representation. 

Now in order to get equations of motion for creation and annihilation oper- 
ators in the OR, let us substitute the decomposition (|^ in Eqs. ||1D|) and use 
relations Then, as an intermediate result, we obtain the following closed 

set of equations of motion which is valid in a general case 

da(p, t) + C/"^(p, t) ai3ip, t) + U^'^i^, t) h^i-p, t) = -itj(p, t) a„(p, i), 

4(p, t) ~ a^(p, t) t/f "(p, t) + 6^(-p, t) utiP, t) = «^(P, t) aj,(p, t), 

6„(-p, t) + at (p, t) "(p, t) - 6^(-p, t) l//"(p, t) = -zc^(p, t) 6„(-p, t), 

6t (-P, t) + V^\y>, t) aMP, t) + ^2"''(P, i) ^M-P, t) = ^^(P, 0&i(-P, 0- (47) 
The spinor construction is introduced here as 

Uf = ut"(p,i)u^(p,i), l^i"'^ =vt"(-p,i)u'3(p,i), 

f/2"'' = ut"(p,i)v'3(-p,i), l/2"^'=vt"(-p,f)v^(-p,f). (48) 

The matrices U2 and Vi describe the transitions between states with positive 
and negative energies and different spins, while the antiunitary matrices Ui 
and V2 show the spin rotations only 

Ul^-Ui, V^^-V2, V^=-U2. (49) 



Equations pT)) are compatible with the canonical commutation relations p6|) . 
Let us write now the u, v - spinors in an explicit form, according to |57j : 

uti(p,i) =yl(p)[cj+,0,p3,P-], 

ut2(p,i) = A(p)[0,W+,p+,-p3], 
vti(-p,t) = A(p) hp3,-P-,^+,0], 
vt2(-p,t) = A(p)[-p+,p3,0,c^+], (50) 

where p± = pi ± ip2, tiJ+ = w + m(t) and A(p) = [2a)aj_|_]~"'^/^. Spin rotation 
matrices (|48p in this representation are equal to zero 

C/i = - 0. (51) 

For the remaining matrices we have J72 = — Vi = t/, where U is the 
hermitian matrix 



mit) 



P3 P- 



(52) 
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Thus, the system of equations of motion ([17|) reduces to the foUowing one: 
aa{p,t) + Uaf:i{p,t)bl{-p,t) = -iu;{p,t)aa{p,t), 
ba{-p,t) - a^^{p,t)Ui3a{p,t) = -iuj{p,t)ba{-p,t). (53) 

3.2 Kinetic equation 

Equations of motion (|53[) do not contain the spin rotation matrices (|51l) and 
they are similar to Eqs. therefore, the KE derivation meets no problem 

now. To be specific, let us introduce the one-particle correlation functions 

9ai3iP,t) = {0\al{p,t)aa{p,t)\0), 

~9a0{p,t) = {OM-p,t)bi{-p,t)\0). (54) 

The differentiation of with respect to time leads to the following matrix 
equations 

g(p,t) = -C/(p,<)G(p,i)-Gt(p,i)C/(p,i), 
~g{p,t) = G(p,t)C/(p,t) + C/(p,i)G^(p,i) , (55) 
where the auxiliary function was introduced 

G„^(p,i) = (0|a)3(p,t)6t (_p,t)|0). (56) 
Together with Eqs. (j55|) . the corresponding equation of motion 

G(p,t) = U{p,t)gipA)-~gipA)U{p,t) + 2tu;{p,t)G{p,t) (57) 

forms a closed set of equations for correlation functions. With the help of Eqs. 
((57|) . one can exclude the auxiliary correlator from the system ([55]) 



r 

<?(p, t) = 2C/(p, t) J dt'[~gip, t')C/(p, t') ~ C/(p, t')gip, t')] cos 20(p, i', t), 

to 

t 

^(p,i) = 2 y dt'[C/(p,i')5(P,i')- 5(P,i')f/(P,0]C/(p,i)cos2(?(p,i',t), (58) 

to 

using zero initial conditions. The subsequent transformation is based on the 
relation 

Tr{U{t)AUit')} = iA(p, t)X{p, t')TrA, (59) 



for an arbitrary second-rank matrix A, which follows from Eq. (|52p . The 
function {p = |p|) 
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plays a role of some analog of Eq. . 

Using isotropy of the considered system, we will limited om'selves to the 
spin-averaged scalar distributions 

f(p,t) = ^Tr g{p,t), f{-p,t) = l-^Trg{p,t). (61) 

Calculating the trace of we get 

t 

f{p,t) - h-P,t) - ^A(p,i) J dt'X{p,t')[l-f{-p,t') 

to 

- f{p,t')]cos2e{p;t,t'). (62) 

In the case of the vacuum initial state, /(— p, t) — /(p, t), we have 

t 

/(p,<) = 2A(p,t) J dt'X{p,t')[l-2f{p,t')]cos[2e{t,t% (63) 

to 

The KE's (|23p and (|63p are similar but differ by statistical factors 1 ± 2/ (the 
Bose enhancement or the Fermi suppression) and by structure of the factors 
pop and (j60p . The corresponding linear equations for the non-Hamiltonian 
dynamical system become 

/ = ^\u, u = A[l - 2/] - 2ujv, V = 2uju. (64) 

This system possesses one first integral of motion (see [55] ) 

{l-2ff + + = 1. (65) 

This relation represents an ellipsoid in the phase space of {f,u,v) variables. 
After exclusion of the function / from Eq. ([64|l . we obtain the system of non- 
linear equations 

ii — \\/ 1 — u'^ — v'^ — 2u}v, 

i) — 2lju. (66) 

It can easily be proved that the KE (f63|) is invariant with respect to time 
inversion. Equations analogous to Eqs. pS)) and ([M]) were obtained in [S3] 
(see also ^5^) for the conformal flat space-time. The KE ([55)1 to the case of 
spinor QED was generalized in JA^Ml- the general case, the spin correlation 
functions ([54| and ([56]) can be decomposed in respect of the Pauli matrices {e.g., 
[58]). 
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3.3 Observables and regularization 

The total particle number density and energy density in the considered case are 
distinguished from the corresponding expressions ([50)1 and (|33p for the scalar 
system by the spin degeneration factor g = 4 (for an equal number of particles 
and antiparticles) 

n{t)^Aj[dp]f{p,t), (67) 
e{t) = (0|Too|0) = 4 J [dpMp, t)/(p, t) , (68) 
where Tqq is the zero component of the energy-momentum tensor 



o[V^7m(5.V')-(5.^)7pV]. (69) 



By definition, the entropy density of the fermion system is equal to 

S{t) = ~4 / [dp] [f In / + (1 - /) In (1 - /)] . (70) 



Finally, the pressure is 

P{t) - i < Tfcfc > . (71) 
Using ((69| and ((40)) . this relation can be reduced to the following form: 

Pit) = ^{e{t)~m{tmx)4,{x))}. (72) 



Here the correlation function is calculated by means of relations (|44)) 
1 If m?{t) 

where the last term takes into account the contribution of the vacuum polar- 
ization, 

Ppoi{t) = -lni{t) [ [dp] — ^ u{p,t). (74) 

This result was obtained under additional conditions for "observable" correla- 
tion functions 

(0|at (p,t)a^(p',t)|0) = (0|6i(p,t)&^(p',t)|0) = S^^pSpp, , (75) 

which is a consequence of space-homogeneity and isotropy (absence of the spin 
moment) of the system. The auxiliary correlation function (j56p is not connected 



16 



with the spin moment and, therefore, it remains ofF-diagonal with respect to 
spin indices. 

The same regularization procedure can be realized here for the calculation 
of divergence integrals, as presented in Sect. l2.3l for the scalar bosons. Equation 
(|25p for the distribution function in the low-density approximation is valid also 
after the replacement A(p, t) A(p, t). Asymptotics of the factor (pO)) is equal 
to A°°(p, t) = m{t)/p. Using Eq. ((25|) and the rules of Sect. 12. 3i one can derive 
the following expressions for counterterms for the case of fermion fields: 

rn{t) 



However, these counterterms can be ignored in computer calculations. 

It is known [26] that the phase density of pairs created in an electric field for 
the whole period of its action is related to long-time asymptotics of solutions 
of some oscillator equations. For the inertial mechanism, analogous derivation 
results in the following relation [52] 

lim /(p t) = gQshr^[(m, - m/)] - coshr^[(cUj - Uf)] ^^^^ 
t-*+oo ' 2 sinh (rTrO^i) sinh (r^tj/) ' 

where T,r — ■nr /2, the indices i, / correspond to the initial and final states. This 
relation is convenient for calculation of observables for long pulses r ^ to, e.g. 
Fig. m where the direct solution of kinetic equation is a very robust numerical 
problem. 



3.4 Numerical results 

Here numerical investigations of the KE's are presented for bosons ([5^. fermions 
and appropriate densities of observable variables are estimated (Sects. 12.31 
and 13. 3p . As an example, two variants of time-dependent masses are consid- 
ered. The first case qualitatively corresponds to a typical meson mass change 
under the phase transition within the NJL model [59j 

TO(t) = (too -TO/)exp[-(t/r)^] -hTO/, i>0, (78) 

with the parameters too (initial mass), to/ (final mass), and r (transition time). 
Another variant suggested in 22 , 23J allows an analytical solution of the Dirac 
equation 

™W = ^^+(^^)tanh(2t/.). (79) 
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Figure 1: Time dependence of the Figure 2: The pair density evolution 
pair density for bosons and fermions of fermions for various initial masses 



at mo — 1. 



(p - proton, e - electron). 
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Figure 3: The residual pair density Figure 4: The residual energy density 
versus the mass ratio. versus the transition time r, mo = 1. 



Numerical results for solution of the KEs are presented in Figs. [T][8]for m{t) 
defined by Eq. ([78|) with r = 10 fm/c and in Figs. IQIlOl for the mass ([79| with 
T = 20 fm/c. Masses are specified in natural units, m/i — 197 MeV/c^. 

At a glimpse the time dependence of quasiparticle density for bosons and 
fermions repeats qualitatively the curve m(t); however, when m{t) — )■ the 
densities go asymptotically to certain finite values (residual density. Fig. 
[U t ^ r), which characterize the real (free) particles (at the active stage 
of the process, m(i) 7^ 0, one may talk about quasiparticles only). The n(t) 
dependence for fermions on the initial mass value is shown in Fig. [2]in the range 
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Figure 5: Entropy density evolution 
of bosons for different values of the re- 
laxation time T, mo = 1. 
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Figure 7: Distribution functions at 
the time i = t = 10 fm/c, jtiq = 1. 
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Figure 6: The distribution function of 
bosons in time-momentum variables, 
Too = 1, t = 10 fm/c. 
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Figure 8: Asymptotic form of distri- 
bution functions at i 3> t. 



from the electron mass to proton one. This dependence is nonmonotonic: With 
increasing too the residual density reaches the maximum at too ~ 10 MeV and 
then begins to decrease. This effect appreciably depends on the variant used 
for the mass change: For the case (f79| it manifests itself much more clearly 
than for the model ([75]) . The residual particle density dependence on the mass 
change is presented in Fig. [3] Qualitative behavior of the energy density 
is similar to that of particle density, showing a smooth asymptotic decrease 
(see Fig. . In Fig. [5] the time dependence of the boson entropy density is 
presented for different values of the relaxation time. Non-monotonic behavior 
of entropy is caused by the fact that the system is opened and treated in the 
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Figure 9: Pressure evolution corre- 
sponding to mass changing ([79|) , toq = 
1, T = 20 fm/c. 



Figure 10: Asymptotic oscillations of 
pressure at i 3> t without 'switching 
off' of vacuum polarization effects. 



non-dissipative approximation. 

Momentum spectra of particles at different stages of the interaction process 
are shown in Figs. [6][8l The maximal number of bosons is created with zero 
momentum, whereas there are no fermions with p = 0. This feature differs 
qualitatively from the case of the Schwinger mechanism of particle creation 

It is important that the formation of appreciably non-monotonic distribu- 
tions with the "fast" mass changing ((78)) . Fig. [3 assists in the development 
of plasma oscillations. For smoother mass changing ([79]) this effect becomes 
much less pronounced. 

The most interesting features are observed in the pressure behavior, Fig. 
[MTUl For both variants of the mass evolution and independently of particle 
statistics, the pressure is negative at the beginning of process, then it changes its 
sign in the reflection point of m(t) and gradually decreases. However, contrary 
to other observables, the pressure has no constant asymptotics and at i ^ r 
looks like almost un-damped oscillations, Fig. [9l It is distinctive for pressure of 
the bosonic quasiparticle system which strongly oscillates around zero. Fig. (TUl 
The reason is that unlike the other considered quantities, the pressure is not 
completely determined by the quasiparticles distribution function /(p,t) , but 
it depends also on the function u(p,i), which describes vacuum polarization 
effects. At the operator language, this means incomplete diagonalization of the 
energy-momentum tensor in the Fock space: Averaged over the initial vacuum, 
its spatial components include the contribution of anomalous correlators like 

(014410). 

Thus, if the process of particle creation stops when the time mass evolution 
is completed {m{t) — > 0), the vacuum polarization effects are not "switched off' 
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simultaneously but continue to influence some observables, e.g., pressure. As a 
consequence, in such non-dissipative nonequilibrium model it is impossible to 
determine unambiguously the equation of state [5U] . 

4 Massive vector bosons 
4.1 The complete QPR 

The simplest version of quantum field theory of neutral massive vector bosons 
is given by the Lagrangian density [52] 

C{x)^-]^d^u,d^'u'' + ]^m^{t)u,u\ (80) 

which corresponds to the equation of motion 

[a^5^ + m2(t)]7i, = (81) 
with the additional "external" constraint 

= 0. (82) 
An alternative way is to proceed from the Wentzel Lagrangian [61] 

C{x) - --^F^'^F^, + ^m^t)u,u'', (83) 

where the strength tensor — d^u^ — d^j^Uy. Here the constraint ([5^ is 
a consequence of dynamical equations. Lagrangians (|80p and ([83]) result in 
different energy- momentum tensors {e.g., see |6H I62j). 

The transition to the QPR is carried out by the standard decomposition of 
free fields and momenta in the discrete momentum space with the replacement 
m — > m(i) in the dispersion law (see Sect. 12. ip . 

u,{x) = -L= ^ ^}—- \a,{v. t) e^^- + a;{p, t) e-^^A, 
V2F ^ V^(p,t) I J 

^^(x) = --L= ^ ^/Ui^{a,ip, t) e^P- - a;(p, t) e^P^} , (84) 

where a^(p,t) are the classical amplitudes. Unlike the scalar case, the consis- 
tent quantization is possible only after including the constraint (|82[) . 
The substitution of the field operators (IM)) into the Hamiltonian 

H = 1 dx (tt^tt'^ -I- Vm^Vu^ + ni^(t)Uf,u'') (85) 
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gives directly the diagonal form in the Fock space 

H = - tj(p, t) a* (p, t) a^(p, t). 



(86) 



However, this quadratic form is not positively defined. To correct it, one needs 
to exclude the /x = component by the additional condition The equa- 

tions of motion for amplitudes a^(p, t) are similar to the scalar case (|15p 



1 



ap(P' t) = 77 A(p, t) a* (-P, t) - iuj{p, t) a^(p, t), 



(87) 



where A(p,t) is defined by Eq. Using Eqs. ([57)) . the condition ([5^ may 

be transformed now to the following relation (i = 1,2,3) 

w(p, t) ao(p, t) = K ai(p, t). (88) 

This equation allows one to exclude the fi — component from the Hamiltonian 
(|86p and quantize other components, which gives 



Pipk 



a|(p,t) afe(p,i). 



(89) 



with the corresponding commutation relation and the vacuum state 



a,(p,i),4(p',i) 



= S,kSpp' a,(p,t^ -c»)|0 >= 0. (90) 

by means 



The next step is the diagonalization of the quadratic form 
of the linear transformation [56] 

ai(p,t) ^Eikak{p,t) = (ei)iai(p,t) + (e2)iQ!2(p, i) + (e3). 



m(t) 



«3(p,i), (91) 



where {ei(p), e2(p), e3(p)} is the local orthogonal basis constructed on the 
vector 63 — p/\p\- These real unit vectors form the triad, 

^ik^jk — ^ki^kj — ^ij 7 ^ik — {^i)k • (92) 

The transformation (|9ip establishes the positively-defined Hamiltonian 



7/ = ^a;(p,i) 



(93) 



The presence of the uj/m factor in the non- unitary matrix E in Eq. ((9T|) leads 
to violation of the unitary equivalence between the a-representation (|89|) and 



22 



a-representations ([M]) . Equations of motion for these new amplitudes follow 
from a combination of Eqs. (jST]) and (|9ip 
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aiip,t) = ^A{p,t)al{-p,t) - iuj{p,t)a^{p,t) + r]ij{p,t)aj{p,t). (94) 



The spin rotation matrix rjij is defined as 



(95) 



with Am ~ —rfi/m + A. This relation shows a particular role of the third 
component. 

Together with the Hamiltonian ((86| . the total momentum operator takes 
also the diagonal form. However, the spin operator 



UfcTTj + TTjUk - UjTTk - TTkUj 



has non-diagonal terms in the spin space in terms of the operator 



al(P:^)aj(P)*) - a{P.t)a]{p,t) 



(96) 



(97) 



where e^-fc is the unit antisymmetric tensor. In particular, the spin projection 
on the p3-axis is 



a\{p,t)a2{p,t) - aj(p,t)ai(p,i) 

+ a2{p,t)a[{p,t) - ai{p,t)a\{p,t) 



(98) 



Thus, this representation can be named the incomplete quasiparticle repre- 
sentation since the spin projection is not fixed. The operator ([M]) can be 
diagonalized by a linear transformation to the new basis [56) 



Ci{p,t) = i?ifcafc(p, i), 



with the unitary matrix 



V2 



1 i 
-i 1 
V2 



(99) 



(100) 



As a result, the new amplitudes Ci{p, t) correspond to creation and annihilation 
operators of vector quasiparticles with the total energy, 3-momentum and spin 
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projection into the chosen direction, 



H{t) = ^ c^(p,i) cl(p,i)cz(p,i)+c.(p,i)cl(p,t) 
p '- 

U{t) = ^p cl(p,i)c,(p,i)+c,(p,i)cl(p,t) 
p '- 

S3{t) = ^ c\{p,t)ci{p,t)-ci{p,t)c\{p,t) 
p '- 

+C2(p,t)4(P,i) - 4(P.^)c2(P,i) 



(101) 
(102) 



(103) 



This c-representation will be referred to as the complete quasiparticle represen- 
tation. The equations of motion for these amplitudes follow from Eqs. (f94 |) .(l99 p 

c,(p, t) = iA(p, i)c|(-p, t) - iuj{p, t)c,(p, t) + ?7y (p, t)cj (p, t), (104) 



where the matrix rjij is fixed by Eq. (|55)) . 

The transition to this representation from the initial a-representation is 
defined by the combination of transformations (PT|) and ([M)) 



c{p,t) = U{p,t)a{p,t) 
with the nonunitary operator (e'^' = (ei ± ie2)/\/2) 



Uip,t)^R-E-\p,t) 



.(+) 



= (-) 
7^31 



7^32 



= (+) 
-3 

J") 



fe33 



(105) 



(106) 



To solve the quantization problem, the equation of motion should be taken 
into account. The commutation relation has the non-canonical form 



[c.,(p,t),c + (p',i)] = Qtk{p,t)Qjk{p,t)Spp', 
where the matrices Qii{p, t) are defined by the equations 

Qtj{v,t) = Vik{p,t)Qkjip,t) 
with the initial conditions 

lim Qtj{p,t) = 6,j. 



(107) 
(108) 
(109) 



So the commutation relation is transformed to the canonical form only in 
the asymptotic limit t — > — oo. Relation (|107p provides the definition of 
positive-energy quasiparticle excitations of vacuum to be treated as some time- 
dependent energy reservoir . 
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4.2 Kinetic equations 

The standard procedure to derive the KE is based on the Heisenberg-type 
equations of motion ([57]) or (|104p . Let us introduce one-particle correlation 
functions of vector bosons in the initial a-representation 

F.kip.t) = (0|a,|(p,i)afe(p,i)|0), (110) 

where the averaging procedure is performed over the in- vacuum state |26j . By 
differentiating the first equation with respect to time, we obtain 

hk{p,t) = \A{p,t) {$,fe(p,i) + $tjp,t)} , (111) 

where the auxiliary correlation function is introduced as 

<^>^k{p,t) = {0\a,{~p,t)ak{p,t)\0). (112) 

Equation of motion for this function can be obtained by analogy with Eq. 
(jllip . We write out the answer in the integral form 

t 

<f.fc(p,t) = ^ J dt'Aip,t')[S,k + 2F,kip,t')]e^''^P-'''''K (113) 

— oo 

The asymptotic condition ^^^^(p, —oo) — (the absence of quasiparticles in the 
initial time) has been used here. The substitution of Eq. pi3p into Eq. (|llip 
leads to the resulting KE 

t 

F,fc(p,i) = ^A{p,t) J dt'A{p,t')[6,k + 2F,k{p,t')]cos[2e{p;t,t')]. (114) 

— oo 

This KE is a natural generalization of the corresponding KE for scalar particles 
(Sec. 

However, there is a number of problems that are specific for the theory of 
massive bosons: The energy is not positively-defined, the spin operator has non- 
diagonal terms in the space of spin states etc. This circumstance hampers the 
physical interpretation of the distribution function (IllOp . In order to overcome 
this difficulty, it is necessary to proceed to the complete QPR where the system 
has well-defined values of energy, spin etc. The simplest way to derive the KE 
in this QPR is based on the application of the transformation (|105p directly to 
the KE (fTT4)) . 

Similarly to the definitions piop . we introduce correlation functions of vec- 
tor particles in the complete QPR 

U{p,t) = (0|ct(p,i)cfe(p,<)|0). (115) 
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They are connected with the primordial correlation functions (|110p by relations 
of the type 

hk{p,t) - Ulip,t)UkmiP.t)Fr,miP,t), (116) 

As a result, the KE (|114p becomes: 
t 

f,k{p,t) = ^A{p,t) J dt'A{p,t')M,kji{p,t,t')[S,i 

— OO 

+ 2f,i{p, t')] cos 20(p; t, t') - A„,(p, t)[5,^hk{p, t) + 43/.3(p, i)], (117) 

where 

w(t') m(i) . . . ' 

H 77; 777TOi3f'fc30j30;3 



(118) 



and (5^4 = - fefe- 

As was expected, distribution functions /q/3(p, t) satisfy the same KE (|114p 
for a = 1,2. The feature of the complete QPR becomes apparent only in tensor 
components of the distribution function fik{p, t) which contains the preferred 
values of the spin index i and (or) fc = 3. Let us select the KE for diagonal 
components of the correlation function pi5p which has the direct physical 
meaning of transversal (i = 1, 2) 

t 

/,(p,i) = iA(p,i) j dt'A(p,i')[l + 2/,(p,t')]cos20(p;i,i') (119) 

and longitudinal components of the distribution function 
/3(p,t) = -2A™(p,i)/3(p,t) + iA(p,t)^ 

X J di'A(p,i')^^[2/3(p,O + Q(P,O]cos20(p;i,O- (120) 



-00 



Here the shorthand notation fa = fi has been introduced for diagonal compo- 
nents of the matrix correlation functions IjllSp and 



mm . . p^ / „ N 

A=^, A^^-A-il^- (121) 



26 



Longitudinal and transversal distribution functions are connected by the rela- 
tion 

/3(p,t) = Q(p,i)/i(p,i), (122) 

where Q{p,t) is the function entering into the commutation relation for the 
longitudinal bosons, 

[c3(p,i),4(p',^)] = QiP,t)Sp 



"pp' J 



Qip,t) = exp 



m{t) uj{to) 



(123) 



with the corresponding initial values m(io) and uj{to). 

Due to relation (fT22l) . it is sufficient to solve only Eq. (fTTOl) . The KE for 
transversal bosons is transformed from the integro-differential form to a set of 
ordinary differential equations similarly to the cases considered above (Sects, 
andf 



fk 



1 ^ 



Uk 



A(l + 2/fe 



2wm,, 



Vk = 2ujUk- 



(124) 



The general initial condition for all diagonal components of the distribution 
function 

lim fk{t) = lim Uk{t) ^ lim Vk{t) = (125) 



t — > — QO t 

satisfies the following requirement 



lim 'm{t) — tuq, or lim rh{t) = 0. 



(126) 



The main operating characteristic of the vacuum creation process is the 
total number density of vector bosons, 

OO 

ntot{t) = j P^dp [2/i(p, t) + h{p, t)] 

p^dph{p,t)[2 + Qip,t)] , (127) 



where isotropy of the system was assumed, p = |p|. As it will be shown in Sec. 
14.31 the integral (|127p is convergent. It essentially differs from results based on 
the Wentzel Lagrangian ([55)1 . In particular, as was shown in [THl [T71 [32] , the 
Wentzel approach leads to infinite particle density. Authors of these papers 
give some arguments in favor of that allowing for quasparticle collisions will 
result in finite physical quantities and bring the estimated energy budget of the 
Universe to agreement with observation data [T71 [53] . 
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We will consider below the case when the time dependence of the vector 
boson mass is defined by conformal evolution of the Universe. The set of 
equations similar to (|124p was obtained first in [HU] within an alternative model 
Lagrangian ([55)1 in the Friedman-Robertson- Walker (FRW) space-time. The 
quantization procedure is quite ordinary in this approach. It would be useful 
to compare the predictions of these two models in detail. 



4.3 EoS for the isotropic case 

The relations for the energy density and pressure can be derived from the 
energy-momentum tensor corresponding to the Lagrangian (|80p 

= -d^Uad^u" - dt.Uad^u" - g^tyC. (128) 

Taking into account the isotropy of the system, we obtain the EoS for the 
massive vector boson gas [TB] 



e{t) = 2 J [dp]Lo{2 + Q)h, 



P[t) 



ls{t) - ^m" J 1^(2 + Q)/i -t- 5P,ac{t), (129) 



where 5Pyac{t) is the contribution to pressure induced by the vacuum polar- 
ization, 



_|/M [i.o(i.g) 



^2 

' ' Ml. (130) 



In order to prove the convergence of the integrals (|129p , (|130p we investigate 
the asymptotic behavior of the solution of the system of equations (|124p . This 
system can be solved exactly in the limit p ^ m for the case a = 1/2 (the 
parameter a is defined in the Sec. 5) 

/ = 7r^4"' ^4(l + 2/)-2p«, v = 2pu. (131) 

AtH P^ 2tH p^ 

The asymptotic solution of Eqs. (|13ip with the initial conditions (|125p is 
,x v{p,t) sm^p{t-to) sm2p{t-to) 

where mo — 'm{to) — rn'Jj^tj^^^ and = l/mi^. The numerical study of Eqs. 
(|124p shows that the basic features of the solutions (|132p for a — 1/2 are 
conserved also for other a > 0. It corresponds to results of Sect. 12.31 
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According to Eq. (|132p , the particle number density and energy density (see 
Eqs. (fT27)) and p29l) . respectively) are convergent, but the vacuum polarization 
contribution to pressure (jl30p is divergent. Moreover, irrelevant fast vacuum 
oscillations of the pressure are observed here. Such a behavior of the pressure 
for a plasma created from vacuum is not specific for the present theory but it is 
inherent in the models where an electron-positron plasma is created in strong 
time-dependent electric fields as investigated in [64]. The standard regular- 
ization procedure of similar integrals with some unknown functions satisfying 
ordinary differential equations is based on the study of asymptotic decompo- 
sitions of these functions in power series of the inverse momentum, (see 
Sect. 12. 3p . In the considered case, such a procedure is not effective because 
the solution (|13ip has quickly oscillating factors (" Zitterbewegung" ), whose 
asymptotic decomposition leads to secular terms. Therefore, for numerical cal- 
culations we regularize the pressure by a momentum cut-off aX p = lOmo and 
separate its stable part by means of the time averaging procedure 



Such "coarse graining" procedure was proposed in [6F in order to exclude 
the "Zitterbewegung" from the description of vacuum particle creation. In 
reality, these fast oscillations are smoothed out due to dissipative processes 
which are not taken into consideration here (the first attempt to derive the 
collision integral for the scalar quasiparticle gas in a strong electric field was 
made in [42j). 

5 Application to conformal cosmology models 

The description of the vacuum creation of particles in time-dependent gravita- 
tional fields of cosmological models goes back to Refs. [45l [66l [SZl [68] and has 
been reviewed, e.g., in monographs [34| , I69| , [70]. The particularity of our work 
consists in the consideration of vacuum generation of particles at conditions 
of the early Universe in the framework of a conformal-invariant cosmological 
model [TOl I17j. Thus, the space-time is assumed to be conformably flat and 
the expansion of the Universe in the Einstein frame (with metric .g^jy) with 
constant masses rh can be replaced by the change of masses in the Jordan 
frame (with metric g^i,) due to the evolution of the cosmological (scalar) dila- 
ton background field [IH [7T] . This mass change is defined by the conformal 
factor r2(a;) of the conformal transformation 



t 




(133) 



(134) 
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Since mass terms generally violate the conformal invariance, a space-time de- 
pendent mass term 

m(x) — m (135) 
il{x) 

has been introduced which formally keeps the conformal invariance of the the- 
ory [33] ■ In the important particular case of the isotropic FRW space-time, 
the conformal factor is equal to the scale factor, il{x) — a{i), and hence 
m{t) = a{t)mobs, where t is "the Einstein time" and mobs is the observable 
present-day mass. Such a dependence was used, e.g., in Ref. j26j for the FRW 
metric. On the other hand, the scaling factor a{t) is defined by the cosmic 
equation of state. For a barotropic fluid, this EoS has the form 

Pph = (7 - = clsph, (136) 

where Pph and Sph are phenomenological pressure and energy density (in con- 
trast to "dynamical" P and e, see Sect. 14.31) . 7 is the barotropic parameter, 
Cs is the sound velocity. The solution of the Friedman equation for such EoS 
leads to the following scaling factor: 

a(f)-p/37^ (137) 

Kinetics of the vacuum creation of massive vector bosons (Sect. 14. 2p was con- 
structed in the flat Jordan frame with the proper conformal time t which is 
necessary to introduce now in Eq. (|137|) . The transition to the conformal time 
is defined by the relation dt — di/a{i). From this relation and Eq. I|137p it 
follows 

37/(37-2) 

(138) 

The substitution of this relation into Eq. (|137p establishes the mass evolution 
law in the terms of the conformal time 

m(t) = (t/iff)" mw, a = (139) 

where tn = [(1 + a)H]^^ is the scaling factor (the age of the Universe), H = 
70 km/s/Mpc is the the present-day Hubble constant and the W-boson mass 
is taken as mw — 80 GeV. Values of the a parameter for some popular EoS 
are: 7 = 2, a = 1/2 (stiff fluid); 7 — 4/3, a — 1 (radiation); 7 = 1, a = 2 
(dust); 7 < 2/3 (quintessence); 7 = 0, a = — 1 (cold matter including baryon 
mass and dark matter) [TH [73| . 

Due to back reactions and the dynamical mass generation during the cosmic 
evolution the detailed mass history remains to be worked out. The central 
question, however, is whether the number density of produced W-bosons could 
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Figure 11: Time evolution of the number density for transversal ni and lon- 
gitudinal n2 vector bosons with the initial condition mo ■ to — 1 for a — 1/2 
(left panel) and the corresponding momentum distribution at the time t ^ to 
(right panel). 




Figure 12: The dependence of the number density of residual vector bosons on 
the initial time to- a — 1/2 (left panel), a — 1/3 (right panel). 

be of the same order as that of the cosmic microwave background (CMB) 
photons, ncMB ~ 465 cm~^. If this question may be answered positively, the 
vacuum pair creation of W-bosons from a time-dependent scalar field (mass 
term) could be suggested as a mechanism for the generation of matter and 
radiation in the early Universe. 

The numerical analysis of Eqs. ()124|) for massive vector bosons is performed 
by the standard Runge-Kutta method on a one-dimensional momentum grid. 
As one can see in Fig. [TTl the creation process ends very quickly and the 
particle density saturates at some finite value. The momentum distribution 
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Figure 13: The time dependence of the energy density and the mean value of 
pressure (|133|) at t :s> to (left panel) as well as pressure (|129p (right panel) with 
the initial condition mo • to = 1 for a = 1/2. 

of particles is formed also very early when m{t) w mo and frozen in such 
form so later on, for times t ^ to, most of particles have very small momenta 
p <C m{t). The spectrum of created bosons is essentially nonequilibrium; hence 
we should continue the analysis of relevant dissipative mechanisms and other 
observable manifestations of a non-equilibrium state (e.g., CMB photons; in 
this connection, see, for example, [Til [75] ). 

The dependence of the corresponding final value of density on the initial 
time is shown in Fig. 1121 The final density ni of the transversal vector bosons 
with spin projection ±1 reaches a maximum when for very early initial times 
we are close to the birth of the Universe. However, in the same limit, the 
density /a of the longitudinal particles with zero spin projection grows beyond 
all bounds. The choice of the EoS changes drastically the number of created 
particles, thus resulting in values which are too small {a = 1/2) or too large 
(a = 1/3) in comparison with the observed CMB photon densities. In order 
to improve this model, we should use an improved EoS, assuming that the 
barotropic parameter 7 characterizing the evolution of particle masses can be 
changed during the time evolution. Such a time-dependence could be induced 
by action of the back-reaction of created particles on the scalar field. Further- 
more, we could use another space-time model, e.g., the Kasner space-time |76) 
instead of the conformally flat de Sitter one. As compared to the earlier work 
|16j . the main achievement of this approach is that there is no divergence in 
the distribution function; thus we do not need to introduce any ambiguous 
regularization procedure. 

As shown in Fig. I13[ the mean pressure remains negative and its magnitude 
becomes negligible in comparison with the energy density (as to the role of 
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Figure 14: The time dependence of 
fermion pair density for a mass evo- 
lution with a ^ 1/2 and niq = 170 
GeV. 



Figure 15: The final density of 
fermions for the conformal time de- 
pendence of mass ()139|) as a function 
of the initial time to for a = 1/2. 



the negative pressure in cosmology see, e.g., [77]). These features can lead 
to violation of the energy dominance condition e + P > that corresponds 
to accelerated expansion of the Universe. Such a kind of models is widely 
discussed (e.g., [78]). For large time moments, the energy density grows but 
the pressure stays very small, P ~ 0. The energy growth under condition P ~ 
results in the conclusion that the massive vector boson-anti-boson gas created 
from the vacuum is cold. It can be seen directly from Eq. (|129p that at large 
times e{t) ~ m{t)ntot{t), because of oj{t) ~ m{t). This EoS of the massive 
vector boson gas (e ^ and p ~ 0) corresponds to dust-like matter [79] , which 
would characterize the evolution of the Universe when the vector boson gas is 
the dominant component of its matter (energy) content. At a qualitative level, 
this conclusion is valid independently of the specific choice of the EoS and, in 
particular, in the case of the dust-like EoS. It would be interesting to obtain a 
formula like (|139p as a result of the solution of the Friedman equation with the 
EoS (|129p - (|130p (such a procedure represents the back reaction problem) and 
investigate self-consistently the production of vector bosons in the Universe. 
Let us remark that vacuum creation of massive vector bosons in the FRW 
metric was first considered in [50] . 

In order to investigate features of the fermion vacuum creation in the confor- 
mal model of the Universe, let us use the corresponding basic KE ([63]) and, as an 
example, choose the system of heavy top quarks with the mass niq = 170 GeV 
(with the change rriw niq in Eq. (|139ll ). The values of the parameter a 
are the same as for the vector boson. As one can see in Fig. [TH the creation 
process is completed very quickly and the particle density saturates at some 
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finite value. 

Finally, Fig. [15] shows the dependence of the residual density on the initial 
time to for m(t) given by Eq. p39p . At the qualitative level, the same pic- 
ture will take place for the vacuum creation of neutralino, which can be the 
main component of dark matter (e.g., [HI])- An analogous problem in early 
cosmology was considered in [IS] . 

6 Systems with metastable vacuum 
6.1 Formulation of the problem 

Here a rather general mechanism of the mass formation as a result of self- 
consistent dynamics of mean-field and quantum fluctuations will be considered. 
The separation of the quasiclassical background field is a common procedure of 
different non-perturbative approaches in QFT [52] [33] [HH [83^ . In the framework 
of this procedure quantum fiuctuations can be described by the perturbation 
theory. 

There is a class of physics problems in which the strong background field 
creates particles which in turn infiuence the background field (the back reaction 
problem). In this respect it is worthy to mention such problems as the decay 
of disoriented chiral condensate ^84j , the resonant decay of CP-odd metastable 
states [4] [85] , the pre-equilibrium QGP evolution [29] [36l [37l [86] , the phase 
transition in systems with the broken symmetry [87] etc. 

The construction of general kinetic theory of such a kind for various poten- 
tials is presented in Sect. 16.21 We will derive the closed system of equations 
for the self-consistent description of the back reaction problem, including the 
KE with nonperturbative source term describing the particle creation in the 
quasi-classical background field and the equation of motion for this background 
field. We use the OR to derive the KE. As an illustrative examples, in Sect. 16.3) 
the one-component scalar theory with and double-well potential are con- 
sidered. In these examples, we study some features of the proposed approach. 
In particular, the problem of the stable vacuum state definition and possibil- 
ity to emerge tachyonic regimes is discussed. As a less trivial example, the 
pseudo-scalar sector of the Witten-De Vecchia-Veneziano model will be consid- 
ered. Similar analysis was carried out in some other models of such kind (e.g., 
[4,i85,i88ii89ji9Qj). In this section, we follow paper [Tj. 
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6.2 The set of basic equations 

Let us consider now the scalar field Lagrangian with a self-interaction potential 

where toq is the bare mass and the potential V[^] is an arbitrary continuous 
function with at least one minimum that is necessary for a correct definition of 
the vacuum state. It is assumed that the field $ may be decomposed into the 
quasiclassical space- homogeneous time-dependent background field (/>o(i) and 
fluctuation part (f>{x) 

^x) ^ Mt) + (141) 

In accordance with the definition of fluctuations, we have (ip) — and ($) = (po, 
where the symbol (. . .) denotes some averaging procedure. The background 
fleld 00 (i) can be treated as quasi-classical one at the condition [91] 

|0o| » l/iStf, (142) 

where St is the characteristic time of the field averaging. 

We consider the case of quite small fluctuations in the vicinity of the back- 
ground field. Therefore, the potential energy expansion in powers of 4>{x) can 
be performed 

= V[(bo] + Ri4> + ^i?2<^' + K[0o, 0], (143) 

where 

and V^[0Oj 0] is a residual term containing the higher order contributions to be 
neglected in the current consideration (non-dissipative approximation). The 
decomposition (I143|) can be finite (for polynomial theories) or infinite. After 
field decomposition ()14H) the equation of motion 

9.9^$ -f -f ^^^^ = (145) 
can be rewritten in the following form: 



where the relation 



9^9'^-m2)</, = Q[0o,0], (146) 



[t) = [0o] = + i?2 [M (147) 
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defines the time-dependent in-medium mass. The term in the r.h.s. of (|146p is 

Q2[0o,0] = i^^0^. (148) 

As a result of averaging of Eq. (|146p . the equation of motion for the background 
field is obtained 

0o + "io'/'o + ^i['/'o]+ < Q2[0o,0] >=0, (149) 

where the time independence of the averaging procedure is taken into account. 

The space- homogeneity assumption implies that the function < Q2[0o, (fl > 
in Eq. ()149|) can depend only on time. As it follows from Eqs. (|148ll and p49p . 
the source term in the r.h.s. of Eq. (I146P is exclusively defined by fluctuations, 

g[(/.O,0] =Q2[0O,'/']- < Q2[0O,0] > . (150) 

On the other hand, in a non-stationary situation the field function (j){x) allows 
the decomposition: 

p 

where (f>^^^{p, t) are the positive and negative frequency solutions of the equa- 
tion of motion 

^(^^\p,t) +u;\p,t),p^^\p,t) ^ -Q[(l,o,^;±p] (152) 

with 

u;2(p,t) =m2(t)-f p2 (153) 
and with the Fourier image Q[(j)Q, 4>; p] of the function Q[(j)o, (A], 



= -^^Q[(l)o,(l>;p]e'' 



(154) 



The function Q[(f>o, (f>;p] gives a nonlinear contribution to Eq. ()152p . We 
suppose that there exists a finite limit lim 0^^^ (p, t) — (j)"^^ (p) in the in- 



t- 



finite past and assume that solutions 4)^ ' (p, t) become asymptotically free 
4>^^\p,t) e^*'^"*, where W-(p) = lim (jj(p,t). The existence of the 

t — > — oo 

last limit is based on the adiabatic hypothesis about switching off the self- 
interaction in Eq. (|147p . 
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After the decompositions (|14ip and (|143p the Hamiltonian density is 

Hm^H[<j)o]+Hi[(j)o,4>]+H2[^o,^]+Vr[(j)o,4'l (155) 

where Hq[(I)q] is the background field Hamiltonian, and he terms Hi and H2 
are the Hamiltonian functions of the first and second order with respect to the 
fluctuation field 

i7[0o] = Ho[M + V[M - + \mlcf>l + V[U, (156) 
i7i[0o,0] = 0o0+(mg0o + i?i['/'o])0, (157) 
772[0o, 0] = + \{yct^f + l^n'^^ (f 58) 

The quasiparticle representation for fiuctuations is constructed by means of the 
decompositions ([S]). In order to obtain equation of motion for operator a(p,t), 
let us write the corresponding action with the Hamiltonian p56p - (|158p (this 
way is alternative to the Hamiltonian approach of Sect. 12. ip 

S[(j}]:^ J d'^x {tt^~ Hi- H2-Vr}. (159) 

If decompositions ([5]) are substituted here, we get 

S[4'] = ^J2 J '^^ {*M(P,<)a(P,i)-a(p,t)at(p,i)] 

K(P'*)«^(-P'*)-«(-P'^)a^(P'^)] (160) 
- Lu{v, i)M(p, i)a(P, t) + a(P, i)a^(P, t)] ~ 2K.p[(/)o, 0]} + Si[4>], 

where 5*1 [0] is the part of the action corresponding to the Hamiltonian (|157p 
and V^[^Oi'/'] is the Fourier image of the residual potential term. Variation 
with respect to a(p, t) and subsequent transition to the occupation number 
representation lead to the Heisenberg-type equations of motion (p 7^ 0) 

a(p, t) = iw^(p, t)at(-p, t) - i[H2 + a(p, i)], (161) 

where 

W^(p,t) = 4^ = ^. (162) 

^ t^(p,t) 2^2 ^ 

In Eq. (|16ip the condensate contribution generated by the action part Si[(t>\ 
is omitted because it corresponds to p = (an appropriate mechanism of the 
condensate state with p = and excitations is absent in the present model). 
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Let us introduce the distribution function of quasiparticles, according to 
Eq. ((18]). Using methods of Sect. [2] and Eq. (|16ip we get the kinetic equation 
in the nondissipative approximation, V^[0o,0] — ^ 0, 



ii = W[l + 2/) - 2ujv, V = 2wit, 



(163) 



which is an analog of p6|l with the replace of A by W . To rewrite Eq. (I149p for 
the background field in the non-dissipative approximation, one has to calculate 
the averaged value < in \ (f)'^{x) \ in >. In the space-homogeneous case one can 
obtain 



< in I (j)^{x) I in >- 
Then Eq. pig)l is reduced to 



1 



dp 



w(p,<) 



[l + 2/(p,<) + u(p,i)], 



4>o + ml(l>o + Ri [0o] 



ldi?2 [ dp 
2 d(/)o J Lu{p,t) 



fip,t) + -u{p,t) 







(164) 



(165) 



(the vacuum term is omitted in integrand). 

The KEs (|163p and (|165p form the closed set of nonlinear equations describ- 
ing the back-reaction problem. In the case of v[^q, $] — 0, this set of equations 
directly follows from nonperturbative dynamics and assumption (jl4ip . For the 
description of the particle creation we will use the particle density ([30)1 as well 
as the background field energy and energy of created quasiparticles Cg 



id 



^<Pl + lrnl^l + V{(t>o), 



(2^ 



oj{p,t)f{p,t). 



(166) 



The conservation of the full system energy can be shown analytically. 

The constructed formalism allows the consideration of certain initial states 
at the time t = to: Initial excitation of the background field 4>o{to) and </'o(^o) 
is given under the additional condition that either /(p, to) = or /(p, to) ^ 0, 
where /(p, io) is some initial plasma distribution. 

6.3 Examples 
6.3.1 ^'^ potential 

The separation of the background field (|14ip in the potential 



A > 



(167) 
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Figure 16: Time evolution for the symmetric potential: (a) the mean field; 
(b) the particle density; (c) the energy density; (d) the momentum spectra of 
particles at time moments t = 40 fm/c and t — bO fm/c. Parameter values 
used are: A = 1, mo ^ mh ^ 197 Mev, 0o(O) = 1.2/to-\ 0o(O) = 0. 



leads to the following decomposition coefficients (|144p : 

V[M = 1^4, Ki'/'o^] = A(0o + </'/4)0^ (168) 

Thus, the time-dependent quasiparticle mass of the fluctuating field (|147p is 
equal to 

m^{t) = ml + 3X4:1. (169) 

If A < and the excitation is strong enough, it is possible that a tachyonic 
mode will arise that corresponds to an unstable state The mass (|169p 

determines the factor (|162p 

»-,p,O.A-^. ,170, 
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The KE (|163p with this factor is correct in the nondissipative approximation 
where the residual potential Vr is neglected. 

Let us write down also the equation of motion for the background field (|165p 
in this approximation: 

0o + M2(t)0o + A03 = (171) 



with the corresponding mass to be equal to 



M\t) =m^ + 3A 



a;(p,i) 



/(P,i) + ^"(P) 



(172) 



i.e. the mass of the condensate excitations is defined by both the distribution 
of quasiparticles and vacuum polarization. 

In numerical calculations we apply zero initial conditions for the distribution 
function and nonzero ones for the background field 0o(*o) = 1-2 ( here and 
below we use the units ft = c = 1). Parameter values are chosen by analogy with 
[87] . where the authors offered an alternative method for describing quantum 
systems under action of the strong background field (so-called Cornwall-Jackiw- 
Tomboulis method [33]). As is seen in FigfTBl at the early evolution stage all 
energy is mainly concentrated in field oscillations. For i < 50 the particle 
number density grows slowly. However, density drastically increases at t ^ 50 
and after this time the quantum energy dominates over classical one. 

The case A < (absolutely unstable potential) is associated with a tachyonic 
regime which is realized for high enough excitations when the initial amplitude 
(^o(io) satisfies the condition m§ + 3A0o(io) < 0. 

6.3.2 Double-Vifell potential 

The potential 

= Ia*"* - A > (173) 

leads to the same equation (|171|) for the background field but with a new mass 
(we put here mo = and > 0) 



The factor (fTB^ equals 



dp 



(174) 



where now the quasiparticle frequency (|153[) includes the time-dependent mass 

m2(t) = -^2 _^3A0o- (176) 
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Figure 17: Time evolution for the bistable potential (|173p : (a) the mean 
field; (b) the particle density; (c) the energy density; (d) the momentum spectra 
of particles at time t = 10 fm/c and t = 40 fm/c. Parameter values are 
mo = 0, Ai = ruh, A = 6, 0o(O) = 0.58/m"i, 0o(O) = 0. 

The vicinity of the central point 4'o{t) = is an instable region. In this 
region the group velocity Vg — du}{k)/dk = k/uj{k) is either super-luminous for 
k > kc, where kc is the root of the equation Lu{k,t) ~ or undefined for k < kc. 
Thus, it is the tachyonic region. 

Let us denote the minimum position of the potential (|173p as — ±^o — 
i/i/y/X and put the new origin of the reference frame in one of these points, 
$ = + We discriminate now the background component (/)o and the 
field ^E", i.e. — (j)o + (p. We will omit the sign indices (±), which identify 
the branch of which the fields 4>q and (p belong to. Only rather a small 
range of excitations in the vicinity of the stable points ^± will be considered 
below, |(/)(a::)|, |(?!>o(t)| <C -\/2^o, where ±-\/2^'o are the roots of the equation 
y[\E'] = 0. Using Eq. (jl45|) and methods described in Sect. 16.21 one can obtain 
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the following system of equations of motion: 

4>o + 2(1^^0 + 3A*±02 + + 3A(*± + 0o)(0') + A(03) = o, 
-[9^9'' + ml{t)]cf, + 3A(*± + <^o)[(0') - 0'] + X[{cp^) - 0'] = 0, (177) 
where 

m|(i) = 2/x2 + 3A0o(<?!'o + 2*±). (178) 

In the nondissipative approximation only linear terms should be kept in Eq. (|177p . 
This leads to the KE p63p with the factor defined by the time-dependent mass 

^(P,<) ^ 3A^(0o_+*±) 
2c^(p,t) 2Lo^p,t) 

The mean values (0^) and {(j)^) are calculated either in the minimal order 
of the perturbation theory (for A <C 1) or in the random phase approximation. 
We use the result (|164p for (0^) and {(f>^) = 0. Thus, we have 



00 + A(/io 



(2 + 3(/)o)«'± + 0a 



+ 3A(*± + 0o) / TT^^ [2/(p, t) + ^'(P, t)] = 0. (180) 

Numerical results for the set of parameters corresponding to [57] are presented 
in Fig. [T71 The stationary regime is achieved faster than in the case of the 
symmetric potential. 

Another formalism for describing the strong field problem in the quan- 
tum field system with the potential (|173p is developed by J. Baacke et al (see 
[87] . and papers cited therein). Here the general kinetic approach is devel- 
oped for arbitrary highly-excited nonequilibrium states in the scalar QFT with 
self-interaction admitting the existence of unstable vacuum states. We restrict 
ourselves to the collisionless (non-dissipative) approximation. However, at- 
tempts to go beyond this approximation have been made 42J. As a particular 
example, 0"* and double-well potentials were investigated. 

It would be of interest to study some other properties of the considered 
model, such features as excitation transitions between states with different 
vacua (in the same space-time point). Apparently, it is possible at high initial 
excitation |0o(O)| > ^/2^'(0) and as a consequence of the tunnelling process 
through the central barrier. The last problem is especially interesting in the 
generalized double-well potential model with non-degenerated vacuum states. 



6.3.3 ?7-meson system 

We will use now the developed technique to consider a more complicated 
quantum field system, namely, the 77-meson sector of the Witten-DiVecchia- 
Veneziano model [2^ which describes the low-energy dynamics of the nonet of 
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pseudoscalar mesons in the large Nc limit of QCD. We drow our attention to 
the singlet state of this model with the following Lagrangian [85] : 



/: = ^a,^a-ry + 5Vcos(|)-^r?^ 



(181) 



where b = ^3/267r, b^ = 92 MeV is the semi-leptonic pion decay constant, 
^2 — i(m^ + 2m|.) ~ 0.171 GeV^, ttItt and rriK are tt- and if-meson masses, 
a = rn^ + mi, - 2m\ - 0.726 GeV^ for zero temperature. The corresponding 
total Hamiltonian density is given by 



H 



1 



1 



1 



-arj 



26^ /i^ sin 



2 Jl 
2b' 



(182) 



where the constant addend was discarded. 

The Hamiltonian density (|182p can be reduced to the form [2] 



1 



Hi, 



2b^n^ 



■ 1 V 
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where 



. The effective potential Hin is constructed in such a way 



that its formal decomposition with respect to the field function ri{x) does not 
contain the corresponding squared contribution to be associated with the mass 
term. However, the redefinition ()183p leads to the absolutely unstable potential 
(Sec. 16.3. ip with the corresponding tachyonic modes 2J. In our opinion, these 
modes have artificial character and can be eliminated by return to the original 
Hamiltonian density (|182p . 

By analogy with Eq. ()14ip . let us select the quasiclassical field r]o{t) and 
quantum fluctuation part (t>{x), 



■q{x) = T]a{t) + (j){x). 



(184) 



Now we can implement the general formalism of Sec. 16.3.21 to the 77-meson 
system with the Hamiltonian density (|182p . The master equation (|165p is now 
given by 
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where w(p,t) is the quasiparticle energy (|153p with the mass (|147p 



m^(t) = a + /x^ cos 



m 



(186) 
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Figure 18: The time evolution for the potential ()183p : The particle density (left 
panel) and the mean field energy (right panel). Parameter values are mo = 0, 
%(0) = 0. 



The KE (|163p is defined via the factor (|16'2p which can be presented as 

4:buj^[p,t) b 

The mass formula (|186l) allows one to fit parameters a, jj? and m?{t) > for 
any amplitude of the quasiclassical field 770 (t), as is seen in Fig. [THl The qual- 
itative behavior of the presented observable is similar to that for the bistable 
potential case. Oscillations of the condensate energy are more pronounced for 
the symmetric potential that is related with higher initial values of 770. 

It was shown that if one vacuum state among a set of vacuum states is fixed, 
it leads to some specific evolution of the system under action of the mean field. 
But there is an open question concerning occupation of other system states due 
to tunnelling process [70j . This problem follows also from our treatment. 



7 Summary 

This review is mainly based on the results obtained by the authors in the last 
years and aims to summarize the information about the kinetic description of 
vacuum particle creation. The latter results for the inertial mechanism with the 
time-dependent particle masses are stipulated at the phenomenological level. 
Three basic quantum field models were considered here: Massive scalar, vector 
and spinor fields. The constructed kinetic theory was applied (Sect. [5]) to a 
conformal cosmology model for investigation of matter created from vacuum 
in an early period of the Universe evolution. In particular, it was shown that 
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the density of the produced vector bosons is sufficient for explanation of the 
present-day density of CMB photons. 

The obtained results can be used for a subsequent study of different aspects 
of matter dynamics created from the vacuum in the early Universe (the equation 
of state, the long wave-length acoustic excitations, the back-reaction problem 
etc.). Only the single mechanism of a mass change was considered where this 
change is induced by the conformal expansion of the Universe whose action is 
switched on for the mass mo(io) at some arbitrary initial time tg. For construc- 
tion of a more elaborated theory, one should eventually take into account the 
inflation mechanism of mass generation acting during an earlier period of the 
Universe evolution [751 [Ml [Ml [SS] ■ It would be of interest to consider the gen- 
eration of particles of different masses and quantum statistics using Eq. (|139p 
which is valid for all particles independently of their inner symmetry. 

Interesting perspectives are opened up also by results of Sect. [H where some 
simplest quantum fleld models are investigated for scalar fields with various self- 
interactions and a corresponding quasiclassical nonstationary field (the problem 
of a phase transition at the restoration of broken symmetry, particle tunnelling 
between states with a different vacua etc.). 

Certainly, the considered examples do not exhaust all variety of systems, 
where the vacuum generation of particles is induced by the inertial mechanism. 
In particular, a large class of experimentally controlled models of meson and 
quark subsystems evolved close to a phase transition has remained beyond our 
scope. In addition, the investigation of two-particle correlations is of interest 
for description of some delicate experimental effects in nuclear reactions [23] . 
The extension of the mean-field approximation to take into account the back- 
reaction and collisions is studied intensively with different methods in some 
scalar models [551 ISZ]- The important problem of relevant initial conditions 
also attracts attention in the context of heavy-ion collisions [97]. We did not 
worry much concerning relativistic invariant form of the developed formalism 
: it is an open question. The basis for its solution is the covariant Hamiltonian 
formalism in relativistic kinetic theory [T31 [311 [32] . 

Finally, some effects discussed here take place also in condensed matter 
physics and can be interpreted in the conceptual framework of the inertial 
mechanism {e.g., [98l [99j and references therein). Such a kind of analogy is 
very promising for an experimental test of the given theory with more accessible 
realizations under condensed matter conditions. 
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